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ABSTRACT 

The initial radial density profiles of dark matter halos are laid down by gravitational collapse in 
hierarchical structure formation scenarios and are subject to further compression as baryons cool and 
settle to the halo centers. We here describe an explicit implementation of the algorithm, originally 
developed by Young, to calculate changes to the density profile as the result of adiabatic infall in a 
spherical halo model. Halos with random motion are more resistant to compression than are those in 
which random motions are neglected, which is a key weakness of the simple method widely employed. 
Young's algorithm results in density profiles in excellent agreement with those from iV-body simulations. 
We show how the algorithm may be applied to determine the original uncompressed halos of real galaxies, 
a step which must be computed with care in order to enable a confrontation with theoretical predictions 
from theories such as ACDM. 

Subject headings: galaxies: formation — galaxies: kinematics and dynamics — galaxies: halos — dark 
matter 



1. INTRODUCTION 

Initial density fluctuations in the early universe are 
widely believed to seed the collapse of dark matter ha- 
los through gravitational instability. Computation of this 
process has allowed the clustering and properties of the 
collapsed halos to be predicted in detail in simulations of 
the collisionless component only (e.g., Navarro, Frenk & 
White 1997; Jing 2000; Bullock et al. 2001; Diemand et 
al. 2004; Navarro et al. 2004). 

Initially, the dark matter and hot gas are well-mixed 
(Spergel et al. 2003), but galaxies are believed to form as 
the baryons cool and settle towards the centers of collapsed 
halos, as originally envisaged by White & Rees (1978), Fall 
& Efstathiou (1980), Gunn (1982), and others. The set- 
tling of gas towards the center causes further compression 
of the halo. This latter process is followed directly in sim- 
ulations that include gas cooling (Gottlober et al. 2002; 
Abadi et al. 2003; Governato et al. 2004), but an approxi- 
mate analytic treatment is needed for many applications. 
The conventional procedure (Blumenthal et al. 1986) is 
widely used to compute, for example, the expected surface 
brightness (e.g., Dalcanton, Spergel & Summers 1997; Mo, 
Mao & White 1998), more-elaborate semi-analytic galaxy 
formation models (e.g., Cole et al. 2000; van den Bosch 
& Dalcanton 2000), predictions for lensing (e.g., Keeton 
2001), predictions from tilted power spectra (e.g., Zenter 
& Bullock 2002), etc. It has been reported that the Blu- 
menthal formula generally over-estimates the actual com- 
pression (Barnes 1987; Sellwood 1999; Gnedin et al. 2004). 

However, Young (1980) showed that adiabatic compres- 
sion of a spherical system can be treated exactly, without 
having to run expensive simulations. His formulation was 
originally to model the growth of a black hole in a spher- 
ical star cluster, but the method can be applied for any 
adiabatic change to a spherical potential. It was first used 



for halo compression by Wilson (2004) and developed in- 
dependently by us. 

We review the standard method for halo compression in 
and describe Young's (1980) method in 321 We give a 
few examples and tests in 21 and illustrate the application 
of the method to data from real galaxies in S|SJ 

2. HALO COMPRESSION 

The usual algorithm for halo compression is easy to im- 
plement, but very crude. It is generally attributed to Blu- 
menthal et al. (1986), although the same idea was devel- 
oped independently by Barnes & White (1984), Ryden & 
Gunn (1987), and others. It assumes that dark matter 
particles conserve only their angular momenta as the halo 
is compressed, which is equivalent to assuming that all 
the halo particles move on circular orbits (see below). In 
a spherical potential, the squared angular momentum of 
a circular orbit is L 2 = rGM(r), and if this quantity is 
conserved as a disk with mass profile M d (r) grows slowly, 
we have 



rMn) = r f [M d (r f ) + (1 - f d )M f (r f )}, 



(1) 



where Mi is the initial total mass (dark plus baryonic) 
profile, (1 — fd)Mf is the desired final dark matter mass 
profile, and r j is the final radius of the mass shell initially 
at radius T{. The quantity f d is the fraction of the ini- 
tial total mass, assumed to be independent of radius, that 
condenses to form the disk. We can substitute for M/(r/) 
by making use of the assumption 



M i (r i )=M f (r s ) 



(2) 



which is sometimes stated that "shells of matter do not 
cross" . We can then find for any desired rf and, through 



© , we can obtain the mass profile of the compressed 
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dark matter halo. For convenience, we denote this the 
Blumenthal algorithm. 

Eq. assumes that the disk mass is taken from the 
halo in equal proportions at all radii. A possible alterna- 
tive assumption is to imagine that the mass of the disk ar- 
rives from an external source, in which case the right-hand 
side reads r/[Md(r/) -I- M/(r/)]. Scaling the compressed 
mass profile by the factor 1/(1 + fd) would conserve total 
mass once more and results in a scaled density profile that 
is not mathematically identical to that resulting from the 
more common procedure, but differences are small every- 
where. Neither hypothesis is likely to be correct, since the 
condensing mass fraction will depend on factors such as 
the cooling rate at each radius, but such refinements are 
unlikely to affect the inner halo density much. 

Flores et al. (1993) report some TV-body tests which in- 
dicated that the Blumenthal algorithm yields a reason- 
able approximation to the compressed halo, and Jesseit 
et al. (2002) reached a similar conclusion. But Barnes 
(1987), Sellwood (1999) and Gnedin et al. (2004) have 
warned that the predicted density profile is more concen- 
trated than found in their TV-body simulations, especially 
in the crucial inner part. Gnedin et al. suggest that the 
above scheme be amended to take account of the radial 
motions of the particles, and recommend the substitution 
M(r)r — > M(r)r in eq. where f is the time-averaged 
radius of each halo particle. Their revised formula gives 
better agreement with their simulations, in which halos are 
not perfectly relaxed, but reduces to the Blumenthal pre- 
scription for equilibrium distribution functions of spherical 
halos. 

A number of authors have reported difficulties when 
trying to apply the Blumenthal method to compute an 
original halo, from the observed compressed form {e.g. 
Weiner, Sellwood & Williams 2001). As long as the adi- 
abatic assumption holds, valid formulae for compression 
can be employed to deduce the decompressed density. 
However, Weiner et al. found that the Blumenthal algo- 
rithm fails for decompression because the radial density 
profile can become multi-valued, which is both physically 
impossible and inconsistent with the central assumption 
that "mass shells do not cross." Valid formulae for adia- 
batic changes must apply irrespective of the direction in 
which the change occurs, and the failure of the Blumenthal 
scheme in this regard is further evidence of its inadequacy. 

3. YOUNG'S METHOD 

As stressed in §3.6 of Binney & Tremaine (1987, here- 
after BT), adiabatic changes conserve all three actions of 
an orbit. In a spherical system, the actions are the an- 
gular momentum and radial action, while the third action 
is identically zero because the plane of each orbit is an 
invariant. Thus adiabatic compression of a spherical halo 
needs to take account of the conservation of radial action, 
as well as of angular momentum. (By ignoring the ra- 
dial action, the Blumenthal algorithm implicitly assumes 
it is zero, and all particles move on circular orbits - as 
is well known.) Young (1980) described an algorithm that 
takes proper account of the effects of random motion in the 
halo, while still assuming spherical symmetry, and simi- 
lar considerations underlie the Fokker-Planck calculations 
of globular cluster evolution (e.g., Cohn 1979). Wilson 
(2004) has already applied the scheme we outline here for 



halo compression. 

The assumption of spherical symmetry is crucial, since 
without it we would need to consider three actions, which 
would make the problem essentially intractable analyti- 
cally - although it could still be followed in TV-body simu- 
lations. Jesseit et al. (2002) and Wilson (2004) tested this 
assumption using ./V-body simulations, and we report an 
additional test in 2]that indeed shows that the spherically 
averaged density profile agrees extremely well at all radii 
with the predicted change from a Young-type code. The 
effectiveness of the spherical assumption results from the 
fact that even quite heavy disks lead to very mild flatten- 
ing of an isotropic halo, which justifies neglect of the third 
action. 

The central idea of Young's method is to apply the con- 
straint that the value of the distribution function (DF) 
expressed as a function of the actions is invariant dur- 
ing adiabatic changes - i.e. fn(J r ,J<f>) — fn(Jr,J<p)- In 
this formula, J r (E,L) is the radial action, and = L 
is the azimuthal action, or total angular momentum per 
unit mass, and the subscripts & n refer respectively to 
the original and new halo profiles. While / expressed as 
a function of the actions does not change, f n (E n , L) does 
change because the relation between the radial action, J r , 
and the specific energy, E n , depends on the potential well. 

Starting from some spherical initial model with density, 
Po(r), potential $oM, an d DF f (Eo,L), we add a second 
mass component, and search iteratively for new functions 
/3 n (r), <E>„(r), and f n (E n ,L). At the n-th iteration, we 
determine a new spherical density profile from 

Pn+ i(r) = 4tt /*" (OC) / LmaX L/ " ( f"' L) dLdE n , (3) 

./<£>„ (r) JO r U 

where the radial speed 

u = {2[E n -^ n (r)]-L 2 /r 2 } 1 / 2 , (4) 

and L max (£ l ) = r{2[E - ^(r)]} 1 / 2 . Since the value of /„ 
to be used in the integrand is determined by the condition 
f n {Jr, J(j>) — fo{Jr, J</>) , we need to know J r (E n , L) in the 
n-th potential well, and then to look up Eq(J Ti L) in the 
initial potential well $>o(r), in order to evaluate fo(Eo, L). 
As these functions are not available in closed form, ex- 
cept for special potentials such as the isochrone (Eggen, 
Lynden-Bell & Sandage 1962), we determine both by in- 
terpolation in 2-D tables of values; the table for J r (E n , L) 
has to be updated at every iteration. 

Since p n +i(r) is spherically symmetric, we have for halo 
compression 

$ n+1 (r) = -G r Mn+ l} r ' ] dr> + $ ext (r). (5) 

J r ' 

Note that <J? ra+ i(r) includes ^cxt, the monopole term only 
arising from the disk mass distribution. It is most con- 
venient to create, and keep updated, 1-D tables for both 
$„(r) and p n {r). 

This new estimate of the gravitational potential is then 
used in eq. to obtain an improved estimate of the halo 
density, and the whole procedure is iterated until the so- 
lution converges. In practice, between ten and twenty it- 
erations are needed for the potential profile to stabilize to 
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a part in 10 . The method conserves mass, the final mass 
of compressed halo differs from the initial by a few parts 
in 10 4 for the grids we typically employ. 

The disk mass can easily be subtracted from the halo 
if desired. All that is needed is to normalize the density 
given by eq. ® by the factor (M/, — Ma)/Mh- Here, M% is 
the mass of the halo to the virial radius. Alternatively, as 
we described above for the Blumenthal algorithm, one can 
re-normalize the combined disk and halo masses to equal 
the original total mass. 

This algorithm is clearly more difficult to program, and 
more time-consuming to run, than the Blumenthal algo- 
rithm. We typically use tables for p n , $„, etc., that have 
100 radial points. The table for J r (E n , L), which has to be 
rebuilt at every iteration, has 100 x 50 values, which when 
fitted with a bicubic spline achieves interpolated estimates 
good to a part in 10 7 or better; constructing this table 
takes the lion share of the computational effort. Never- 
theless, the whole iterative scheme converges within a few 
cpu minutes. 

Young (1980) shows that a given mass profile is less 
compressed when an isotropic DF is assumed than when 
all particles have circular orbits, which is the implicit as- 
sumption of the Blumenthal algorithm. Thus allowance for 
radial motion makes a dark matter halo somewhat more 
resistant to compression. 

It should also be noted that even if the initial DF were 
isotropic, the compressed DF will not be. Since the DF 
of the compressed halo is computed anyway, Young's al- 
gorithm yields direct estimates of the radial anisotropy 
expected. 

The method can, in principle, be adapted for halo de- 
compression. Provided an equilibrium DF for the com- 
pressed halo can be constructed, the algorithm guarantees 
that the decompressed halo must have a positive, single- 
valued, density everywhere. The only significant barrier 
to extending the method to this application is the con- 
struction from observational data of suitable DFs for com- 
pressed halos containing disks. 

Young's method is not restricted to any particular ra- 
dial mass profiles for either the disk or the halo; the added 
external mass can have any spherically averaged mass pro- 
file. Neither is the method restricted to isotropic DFs; any 
f(E,L) that yields the initial spherical mass profile may 
be used. 

4. EXAMPLES AND TESTS 

Figure ^ presents results from an A-body simulation 
designed to test the adequacy of the spherical approxima- 
tion for the disk potential. This test uses the NFW density 
profile (Navarro, Frenk & White 1997) 



Pn(t) 



r(r + r s ) 2 



(6) 



with p s setting the density scale at the break radius, r s . 
We define a mass unit as Mm — Anp s r^ sl which is equal to 

The radial mass 



the mass enclosed within r ~ 5.3054r s 
profile of an exponential disk is 




m{R) = M d 



1 



1 
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Rd 
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Rd 



(7) 
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Fig. 1. — Test of the spherical approximation for a disk grown in 
an NFW halo. The jagged lines show the density profile estimated 
from the particles before and after a thin disk is grown. The dotted 
and smooth solid lines show respectively the density of the ana- 
lytic NFW halo and the compressed profile computed from Young's 
algorithm. The dashed curve shows the compression predicted by 
the Blumenthal algorithm, which is clearly inconsistent with the 
Af-body results. 



where Md is the total mass of the disk and Rd is the disk 
scale length. We characterize a simple disk-halo model by 
two parameters: the mass ratio m d = Md/Misi and length 
ratio £ d = R d /r s . 

The isotropic DF for an NFW halo is given by Edding- 
ton's inversion formula (BT, §4.4), and is positive every- 
where. This equation has to be solved numerically for the 
NFW potential, which requires a sophisticated quadrature 
method. It is convenient to impose an upper energy cut- 
off, so that the halo density declines to zero at a large 
radius, but the boundary term should not be included. 
Alternatively, one can use Widrow's (2000) fitted approx- 
imation for the DF. 

The halo is represented by 1 million particles selected 
smoothly from the DF as described in Appendix B of De- 
battista & Sellwood (2000); in order to reduce shot noise 
when estimating the density in the inner halo, the halo 
particle masses were scaled as y/~L, where L is the total 
angular momentum of the particle. The disk in the simu- 
lation is represented by 100 000 equal mass particles that 
are added at a constant rate in the z — plane and held 
fixed in their initial positions in order that the disk mass 
profile remains exactly exponential. The halo particles 
move in response to forces from both the disk and halo, 
which are computed on a spherical grid of 500 radial shells 
(McGlynn 1984; Sellwood 2003, Appendix A) and include 
all multipole terms < I < 8. The disk is grown over the 
period < t < 250(R d '/GMd) and the final density shown 
in Fig. n is measured after twice this period. The time 
step was 0.005(i^/GM d ). 
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Fig. 2. — Comparison between halo compression using the Blumenthal algorithm on the left, and Young's formula on the right, showing 
that halos are less easily compressed when velocity dispersion is taken into account. In both panels, the solid curve shows the total circular 
speed, the dashed curve is the disk contribution, the dotted curve is the uncompressed NFW halo, and the dot-dash curve is the compressed 
halo contribution. Note that = 0.05 in this case. 



The dotted curve in Fig. ^ shows the density function 
(eq. HjJ) that is closely followed by the initial density profile 
of the 1M particles (jagged line). The smooth solid line 
shows the predicted density profile from Young's algorithm 
after an exponential disk with = 0.1 and = 0.05 is 
grown slowly. The jagged solid line that follows this curve 
closely is the (spherically averaged) density profile mea- 
sured from the particles when the asphericity of the thin 
disk and halo are taken into account. Evidently, the spher- 
ical approximation is quite adequate in Young's algorithm. 

The dashed curve in Fig.^shows the prediction from the 
Blumenthal algorithm, which is clearly inconsistent with 
the A^-body results. This discrepancy confirms the conclu- 
sion reached previously by Barnes (1987), Sellwood (1999) 
and Gnedin et al. (2004) that the Blumenthal algorithm 
over-estimates halo compression in the inner parts. 

A more conventional indication of the difference between 
the predictions by Young's and Blumenthal's methods is 
illustrated in Figurc|21for the same model shown in Fig.^ 
The compressed halo density (dot-dashed curve) predicted 
by the Blumenthal algorithm is shown on the left, while 
the right-hand panel shows the prediction from Young's 
procedure. Note that the greatest differences occur in the 
inner region where the disk contributes, which is generally 
the region most strongly constrained by observed rotation 
curves. 

Figure shows the effect of changing the velocity dis- 
tribution in a Plummer halo model. In all cases shown, 
the Plummer sphere is compressed adiabatically by addi- 
tion of an exponential disk having one tenth the mass of 
the halo, and disk scale length one tenth of the Plummer 
core radius. Dejonghe (1987) provides a family of DFs for 
this model with a single parameter, q, that determines the 
shape of the velocity distribution. The solid curve shows 
the resulting rotation curve when the halo has an isotropic 
DF (q — 0). The dot-dashed curve shows that a DF that 



is maximally radially biased (q — 2) is compressed some- 
what less, while the dotted curve shows the result when 
a strongly azimuthally biased DF (q = —15) is assumed. 
The dashed curve shows the result when the Blumenthal 
algorithm is employed, again showing that its assumption 
of extreme azimuthal bias in the orbits of halo particles 
allows too strong a compression of the halo. 

The compressed density resulting from Young's and Blu- 
menthal's methods generally shows larger differences for 
models with cores, such as King and Plummer models, 
than for cusped density profiles. However, even for the 
NFW profile (Fig. |3J) , the differences remain quite signifi- 
cant. 

5. APPLICATION OF YOUNG'S METHOD TO DATA 

In this section, we use Young's algorithm to deduce the 
initial halos of five galaxies from their observed rotation 
curves. We present these results here purely to illustrate 
an application to real data, and will report a more system- 
atic study in a later paper. 

We selected galaxies with extended HI rotation curves 
from the compilation of Sanders & McGaugh (2002) so 
that the shape of the rotation curve is measured well into 
the halo. In order to explore a range of galaxy properties, 
we selected five galaxies that span a large range in both 
luminosity and surface brightness. The five objects are 
the bright galaxy NGC 2903, the intermediate luminosity 
high and low surface brightness galaxies NGC 2403 and 
UGC 128, and the high and low surface brightness dwarf 
galaxies NGC 55 and NGC 1560. 

We stress that these galaxies were selected to be purely 
illustrative of our technique over a wide range of galaxy 
properties, and do not constitute a representative sample. 
They do, however, provide some interesting insights. 

We assume an exponential disk, which is adequate for all 
five galaxies in this very preliminary study, and an NFW 
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Table 1 

Disk, Compressed, and Primordial Halo Parameters 



Galaxy 


Vf 




y c 


Q 


m d 


(■A 


r b 


c 


V200 a 


Refs. 


NGC 2903 


185 


2.0 


2.2 


0.6 


0.10 


0.12 


16.7 


9.2 


115 


1,2 


NGC 2403 


134 


2.1 


0.6 


0.4 


0.02 


0.15 


14.0 


10.1 


106 


1,2,3 


UGC 128 


131 


10 d 


1.1 


1.0 


0.02 


0.19 


52.6 


3.2 


126 


4,5 


NGC 55 


86 


2.7 d 


0.2 


1.0 


0.01 


0.10 


27.0 


4.6 


93 


6 


NGC 1560 


72 


2.6 d 


1.0 


1.0 


0.01 


0.12 


21.7 


4.9 


80 


7 



a In kins" 1 . 
b In kpc. 

c In solar units, Mq/Lq, in the B-band. 

d This baryonic scale length is rather larger than that of the stars alone. 



References. — 1. Begeman (1987). 2. Wevers, van der Kruit & Allen 
(1986). 3. Blais-Ouellette et al. (2004). 4. van der Hulst et al. (1993). 5. 
de Blok, van der Hulst & Bothun (1995). 6. Puche, Carignan, Wainscoat 
(1991). 7. Broeils (1992). 




Fig. 3. — The effect of velocity anisotropy in a Plummer halo. 
The dashed curve shows the circular speed resulting from compres- 
sion using the Blumenthal algorithm while the other curves show 
the same quantity computed by Young's method for different initial 
DFs. The solid curve is for an isotropic model, the dot-dashed curve 
for a model with maximum radial bias, and the dotted curve shows 
a strongly azimuthally biased model. Note the suppressed zero on 
the velocity axis — all curves continue approximately linearly to the 
origin. 



halo with an isotropic DF. We have to choose values for 
our two parameters: the length and mass ratios, id and 
rrid, of the disk and halo. The model halo is dimension- 
less, so the physical parameters of the primordial halo can 
be inferred only after scaling to fit a particular galaxy. 
The conversion to (c, V200) depends weakly on the Hubble 



constant; we adopt H = 75kms _1 Mpc -1 for consistency 
with the galaxy data. 

We first constructed a grid of compressed NFW mod- 
els with a range of plausible rrid and id and examined the 
final rotation curves from this grid to judge which combi- 
nations of disk parameters might make tolerable approx- 
imations to the data. Crudely speaking, m<j < 0.1 and 
id < 0.2 in NFW halos yield plausible- looking rotation 
curves. Larger rrid compresses the halo too much, leading 
to rotation curves that have a pronounced peak and fall 
too far before leveling off. Larger id causes the disk to 
extend too far out into the halo where the declining por- 
tion of the rotation curve is too obvious compared to data. 
This is only an approximate statement: some individual 
galaxies might be described by more extreme parameters, 
and different halo types could well require very different 
combinations of nid and l&. 

A trend in our grid of models mimics one that is well- 
established by the data. The rotation curves rise quickly 
and fall gradually for larger values of rrid, qualitatively as 
observed for luminous spirals. The disk-driven peak de- 
clines as we employ progressively lower rrid, qualitatively 
giving the gradual, and continuing, rise characteristic of 
dwarf galaxy rotation curves. Varying id at fixed ad- 
justs the amplitude of the total rotation curve relative to 
that of the baryons. 

We used this grid to estimate plausible parameters for a 
first approximation to the data for each case, and then re- 
fine the model iteratively. At this juncture, goodness of fit 
is judged by eye. Our models match the data reasonably 
well, though are certainly not perfect. The shortcomings 
of the fits can in some cases be attributed to deviations 
from our idealized exponential disk model, which we as- 
sume here for simplicity. Indeed, we expect to improve 
the fits by using the exact disk mass profile, but leave this 
refinement to future work, as the dominant uncertainty is 
not in the precise shape of the baryonic mass distribution 



6 



but in its amplitude, which depends on the stellar mass-to- 
light ratio T*. For a specific choice of Y», matching both 
the total and baryonic rotation curves simultaneously re- 
stricts the plausible values of and fairly well, though 
we make no claim that our fits are unique. 

We must adopt some estimate for the baryonic mass 
in each galaxy, which is straightforward for the gas, but 
fraught for stars. McGaugh (2004) discusses the T* for all 
the galaxies considered here, preferring a prescription that 
minimizes the scatter both in the baryonic Tully-Fisher 
relation (McGaugh et al. 2000; McGaugh 2005) and in 
the mass discrepancy-acceleration relation. He defines the 
ratio of the adopted Y„ to this optimal one is given by the 
parameter Q — T*/T opt , which encapsulates the goodness 
of the mass-to-light ratio in this respect: the further Q is 
removed from unity, the larger the scatter induced in these 
otherwise tight correlations. 

The basic galaxy data and results are given in Table 
^ The first column gives each galaxy's name. The second 
column gives the observed flat rotation velocity. The third 
column gives the physical scale length of the exponential 
disk chosen to approximate each galaxy's baryon profile; 
since it is the sum of stars and gas that matters, the scale 
length has been stretched where necessary to approximate 
the total. Column 4 gives our final adopted mass-to-light 
ratio (in the £?-band) for the stars. Where possible, 
this has been held fixed at the optimal value described by 
McGaugh (2004), but in some cases we had to reduce Y* 
in order to obtain a fit, and we give the deviation from the 
optimal value in column 5. This deviation is inevitably 
in the sense that the disk mass must be reduced in order 
to accommodate the cusp of the NFW halo. Columns 6 
and 7 give the halo parameters rrid and id of our adopted 
fit. The scale length of the halo r s corresponding to the 
baryonic scale length Rd and fitted Id is given in column 
8. The NFW concentration and V200 parameters of the 
primordial halo are given in columns 9 and 10. References 
to the original sources of the data used are given in column 
11. 

The choice of stellar mass-to-light ratio has a pro- 
nounced effect on the mass model. There is a tremendous 
range of possibilities between maximum and minimum disk 
(e.g, Fig. 0J. Attempts to grow a maximum disk in an 
NFW halo inevitably lead to rotation curves which peak 
too high and fall too far before flattening out. For this rea- 
son, maximum disks and halos with central cusps appear 
to be mutually exclusive. 

Nevertheless, a heavy disk is important for falling rota- 
tion curves. The disk fraction md modulates the shapes 
of rotation curves. Systematic variation of nid with rota- 
tion velocity appears to be required in order to match the 
observed trend in rotation curve shapes from falling for 
bright galaxies to rising for dwarfs (e.g., Persic & Salucci 
1991). This variation can be rapid, with md ~ 0.1 for 
NGC 2903 falling to m d « 0.02 for NGC 2403 (Fig.®. 

Systematic variation of md has important ramifications 
for the Tully-Fisher relation. It implies that the observed 
range of baryonic disks are all found in a relatively narrow 
range of halo masses. Such a finding changes the slope 
of the Tully-Fisher relation from that nominally predicted 
for halos only. 

UGC 128 (Fig. EJ> is comparable to NGC 2403 in lu- 
minosity and rotation velocity. Despite the low surface 
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Fig. 4. — The rotation curve and mass model of the luminous 
high surface brightness galaxy NGC 2903. The rotation curve data 
are plotted as circles with error bars. The contribution to the rota- 
tion by the baryonic component (stars plus gas) is denoted by the 
square symbols. Large squares are for the modeled mass-to-light ra- 
tio specified in Tabled Also shown as small squares are the limiting 
cases of minimum (gas only with T* = 0) and maximum disk. The 
dashed line shows the adiabatically formed exponential disk used 
to approximate the observed baryon distribution. The solid line is 
the total rotation due to disk plus compressed halo. The primor- 
dial NFW halo (with parameters given in Table HI is shown by the 
dotted line, and the compressed halo by the dash-dotted line. 
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Fig. 5. — The intermediate luminosity, high surface brightness 
galaxy NGC 2043. All symbols and lines are the same as for Fig. ||] 
In addition to the extended HI rotation curve of Begeman (1987), 
the Fabry-Perot data of Blais-Ouellette et al. (2004) are shown as 
the small solid points in the rising part of the rotation curve. The 
extended HI in this galaxy makes it difficult to approximate the 
baryons with a single exponential disk. 

brightness of this galaxy, the halo is perceptibly com- 
pressed: the disk is not completely negligible. The disk 
fraction and £4 are very similar to those of NGC 2403. This 
is hardly surprising, since the rotation curves of the two 
galaxies are indistinguishable when normalized by scale 
length: V(r) is very different, but V(r/Rd) is much the 
same (de Blok & McGaugh 1996). Consequently, the den- 
sity of the halo of UGC 128 must be lower than that of 
NGC 2403, simply because r s is larger at comparable mass. 
This is reflected in a low concentration (c w 3). Though 
the fit with a compressed NFW halo is good, primordial 
halo concentrations this low never occur in any plausible 
cosmology (Navarro et al. 1997; McGaugh et al. 2003). 

NGC 55 (Fig. 0) and NGC 1560 (Fig. ®) are nearby 
dwarf galaxies. The gas is not negligible in these galaxies. 
The baryonic rotation curve is reasonably well approxi- 
mated by a single exponential with a scale length stretched 
relative to that of the stars alone (roughly by a factor of 
two) to represent both stars and gas. 
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Fig. 6. — The intermediate luminosity, low surface brightness 
galaxy UGC 128. All symbols and lines are the same as for Fig. [I] 
This galaxy is well matched to NGC 2403 in terms of luminosity 
and circular velocity: these two objects occupy the same location in 
the Tully- Fisher diagram (de Blok & McGaugh 1996). An approx- 
imation to the baryonic component is obtained by stretching the 
exponential disk scale length to 10 kpc, somewhat longer than that 
of the stars, which varies from ~ 6 to 9 kpc, depending on bandpass 
(de Blok et al. 1995). 
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Fig. 7. — NGC 55 is a high surface brightness dwarf galaxy. All 
symbols and lines are the same as for Fig. |1] The baryonic rotation 
curve is approximated by an exponential disk of scale length 2.7 kpc 
(vs. 1.6 kpc for the stars alone). 
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Fig. 8. — NGC 1560 is a low surface brightness dwarf galaxy. All 
symbols and lines are the same as for Fig. [I] The baryonic rotation 
curve is approximated by an exponential disk of scale length 2.6 kpc 
(vs. 1.3 kpc for the stars alone). 

In order to produce something like the observed, grad- 
ually rising rotation curve, it is necessary to consider very 
small disk fractions. The cases illustrated with m<i = 0.01 
provide a good fit to the baryonic components, and a rea- 
sonable match to the outer parts of the rotation curves. 



The fit to the inner few kpc of the rotation curves is not 
good. The problem here stems from the assumed form of 
the primordial halo. The NFW halo is too centrally con- 
centrated to begin with, a well known problem for dwarf 
galaxies (Cote, Carignan & Freeman 2000; de Blok et al. 
2001; Swaters et al. 2003; de Blok, Bosma & McGaugh 
2003). 

6. CONCLUSIONS 

We have used the algorithm described by Young (1980), 
which takes proper account of radial motion, to compute 
the density of a dark matter halo that is compressed by the 
growth of a disk in its center. Young developed the method 
for compression of a spherical star cluster by growth of 
a central black hole; the algorithm was applied indepen- 
dently for halo compression by Wilson (2004). While 
Young's algorithm is more complicated to program, and 
more time consuming to run, than the simple and popular 
algorithm (Blumcnthal et al. 1986), it correctly predicts 
rather less halo compression, especially in the inner region 
that is most accessible to observations. 

The predicted compressed density profiles agree almost 
perfectly with those found in high-quality iV-body simu- 
lations (Fig. ^) , and are significantly less compressed than 
predicted by the Blumenthal algorithm. The almost per- 
fect agreement with ./V-body simulations confirms also that 
the spherical approximation for the disk mass profile is 
entirely adequate, as previously found by Barnes (1987), 
Jesseit et al. (2002), and Wilson (2004). 

The reason that halos are less easily compressed than 
the Blumcnthal algorithm would predict is that this crude 
method neglects the extra pressure caused by radial mo- 
tions. We present tests (Fig. to show that the com- 
pressed density is lowered as the importance of radial mo- 
tion is increased. Since dark matter halos are formed in 
a hierarchical collapse, the outer halos are likely to have 
radially biased velocity distributions, although the inner 
halos are expected to be close to isotropic. 

The algorithm can be applied to all spherical density 
profiles having valid distribution functions, and any ar- 
bitrary mass profile for the disk. The key assumptions, 
that the halo be spherical and that the disk mass was as- 
sembled adiabatically, are no more restrictive than for the 
usual Blumcnthal algorithm. It is expected that hierarchi- 
cal galaxy formation leads to aspherical dark matter halos, 
and compression and shape changes can be computed in 
these more realistic circumstances only by iV-body simu- 
lation (e.g. Kazantzidis et al. 2005). 

The adiabatic assumption is likely to hold if the disk re- 
mains closely axisymmetric, as the cooling and settling of 
gas is likely to occur over many crossing times in the inner 
halo. However, non-adiabatic resonant interactions will 
develop between rotating, non-axisymmetric structures in 
the disk and the orbits of halo particles (Tremaine & Wein- 
berg 1984). It is unclear at present whether such inter- 
actions can cause significant changes to the halo density 
profile (Weinberg & Katz 2002; Sellwood 2003). 

As proof of principle only, we have derived the parame- 
ters of initial, uncompressed NFW halos for a few well- 
studied galaxies, and plan to apply the technique to a 
larger, and more representative, sample of galaxies in a 
later paper. The initial results are promising in many re- 
spects. We found it easy to obtain reasonable fits, with 



8 



variations in disk fraction naturally explaining the vari- 
ation in shapes of rotation curves. However, it does not 
appear likely that an improved treatment of halo compres- 
sion will provide a resolution of the cusp/core controversy, 
or of the absolute halo density of low surface brightness 
galaxies, which remains a problem for ACDM. 
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Oleg Gnedin, provided a helpful report. This work was 
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JAS and AST-0206078 and NAG513108 to SSM. 
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